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Abstract: 

While the hard phase of the strong interaction is well described by perturbative QCD, the 
soft hadronization phase is less understood. Benefiting from the high statistics from e + e~ 
experiments at the Z° resonance, it is possible to impose strong two-jet cuts on the data 
without loosing the statistical significance. In these events perturbative activity is suppressed 
and hadronization effects can be more prominent. We show that after proper event cuts a 
set of observables are sensitive to differences in the hadronization models. These observables 
can thus be important tools for a more detailed study of the hadronization mechanism. 
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1 Introduction 



High energy reactions like e + e~ — > hadrons are usually described in terms of two phases, 
an initial hard perturbative phase formulated as a parton cascade followed by a soft non- 
perturbative hadronization phase, for which phenomenological models are needed. The first 
phase can be calculated from perturbative QCD, but as exact perturbative results in most 
cases are available only to second order approximate schemes have to be used. In the large 
N c limit planar diagrams dominate |l|], and for e + e~-annihilation, descriptions in terms of 
parton cascades have been quite successful. 

Global features like shape variables (thrust, sphericity, oblateness, major, minor etc.) 
and inclusive particle spectra in e.g. xp and p±, reflect mainly the parton distribution in 
the initial perturbative phase. These observables can be well described assuming that the 
hadron distribution closely follows the distribution of partons (local parton hadron duality, 
LPHD provided a locally invariant cut-off is applied to the cascade. One example is 
the observed peak position of the distribution in the variable £ = hx(l/xp), which for most 
hadrons agrees with parton level results, using a virtuality cut-off Qq which is larger for 
heavier hadrons. Also in some more phenomenological models, e.g. in the cluster fragmen- 
tation model || and the UCLA model [[| the hadronization properties are to a large extent 
determined by the hadron masses. 

A more detailed description of the final hadronic state, including e.g. flavour, the baryon 
to meson ratio, spin and polarization as well as correlations (including Bose-Einstein corre- 
lations), depends also on the non-perturbative properties of QCD, including the structure 
of the vacuum condensate. These features cannot be obtained from perturbative calcula- 
tions. Although some of them may in the future be accessible e.g. from lattice calculations, 
at present it is unavoidable that a model attempting to give a detailed description of the 
hadronization process must contain a significant number of phenomenological parameters. 

To gain insight in the properties of QCD in the soft region and the vacuum structure it is 
essential to isolate the hadronization process from the perturbative cascade. This is however 
not an easy problem. At the end of the cascade the running a s becomes relatively large. The 
result depends on a necessary infrared cut-off and nonplanar diagrams and interference effects 
may be important. This implies that this phase is more uncertain and therefore difficult to 
isolate. It is also possible that the separation between the perturbative and non-perturbative 
regions cannot be well defined. As discussed in ref || a change in the perturbative cut-off can 
to a large extent be compensated by a modification of the hadronization parameters. The 
transition region between the two domains may also show interesting coherence effects [0. 

In this paper we will study ways to isolate high energy systems with little perturbative 
activity, in which the non-perturbative features ought to be more prominent. With the 
very high statistics available from the LEP1 experiments, it is possible to impose event 
cuts which exclude a significant amount of perturbative activity, and still have remaining 
events numerous enough for detailed studies of hadronization. It is possible to calculate 
purity and efficiency measures for different event shape cuts. Our result is however that a 
simultaneous purity and efficiency close to 1 cannot be obtained; there is no cut that can 
be called optimal. Instead different cuts are suitable for a study of different features of the 
hadronization process. 

The outline of this paper is as follows: In section Q we study a set of different event 
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cuts. In section [3] we study some observables which have been proposed for a separation 
of different hadronization mechanisms. In section f| we discuss a few models with different 
properties and in section [| their predictions for the considered observables. Conclusions are 
given in section ||. 

Much of the discussion in this paper will concern different aspects of transverse momenta. 
These will be denoted k± for partons and p± for hadrons. We will also discuss a sum of hadron 
p±, which we will call Q±. 



2 Cuts 

We will here discuss a set of event shape cuts suitable to extract events with little perturbative 
activity. The performance of the cuts are investigated using Monte Carlo simulations, where 
the underlying parton state is known. 

For simulated events we will here with "little perturbative activity" mean events where 
all gluon emissions occur at a k± below some kj_ . (In the analysis presented below we will 
take kj_ = 2 GeV.) The reason for choosing this definition, rather than e.g. a cut on the 
sum of E± for all gluons, is that one hard gluon gives rise to a distinct three jet structure, 
which introduces long range correlations in azimuthal angle. Such a gluon also increases the 
phase space for softer ones, which enhances its effects on the final state. 

Through further emissions, the energy of one gluon is distributed among its emission 
products. As this rather enhances than reduces the importance of the perturbative phase, 
we require k± to be smaller than k±o for each emission, and not for the final partons at the 
end of the cascade. 

We examine the performance of the cuts in terms of purity and efficiency Purity is 
defined as the rate of events where the highest k± for a gluon emission, fcj_ max , is below k± . 
If the purity is low, perturbative gluons may shadow many hadronization effects. Efficiency 
is defined as the acceptance rate among the desired events with k± max <k± . If the efficiency 
is low, the sample might be biased, and the interpretation of the results is more difficult. 
Low efficiency also implies low statistics. 

The cuts are applied to hadronic final states from Z° decays, in events generated by 
the Ariadne and Jetset MC |7|, §. Ariadne is an implementation of the colour dipole 
formalism for a QCD cascade M. This cascade is convenient to use in the present analysis, 
as it is ordered in k±, which implies that k± maiX of the cascade is easily extracted as the k± 
of the first gluon emission. Jetset is a MC for the Lund string fragmentation model |ITJ . 
which in general gives a good description of e + e~ data. 

The considered two-jet measures are 

• 1 — T (T=Thrust=max WT ^j^p )- 
. M (M=Major=maxn M ±n T ^jfp)- 

• &±cius= transverse momentum of three-jet configuration found using fc_|_-based cluster 
method (we have chosen the Diclus algorithm |TT[| ). 

• Sphericity S, which measures the sum of p\ w.r.t. the axis minimizing this sum. 



2 



• "Linearized sphericity" L, which is similar to sphericity, but linear in particle momenta. 

• Multiplicity 

Sphericity S is determined by the largest eigenvalue A max to the sphericity tensor 

_ Ei Pj Pj /j\ 

Ei |pd 2 ' 

through the definition S = |(1 — A max ). Thus S ~ for a two-jet event and 5 ~ 1 for an 
isotropic event. Sphericity is quadratic in particle momenta, which implies that the value 
is changed if a particle is split in two collinear ones. A linear correspondence to sphericity, 
suitable for perturbative calculations |12j, is obtained by a modified tensor 

j^af3 — YliP'iPi /|Pjj 

The quantity L = — A^ x ) thus measures the summed Pj_/|p| w.r.t. the event axis. In this 
paper we will use the terminology suggested by Sjostrand, and call this measure linearized 
sphericity. 

The event shape cuts in thrust, sphericity and jet-fcj. are widely used to obtain two-jet 
events. Apart from these quantities, we have chosen to investigate two measures which can 
be expected to be closely related to the maximal kj_ of the parton cascade. One is major, 
which is given by the maximal summed |p| component orthogonal to the thrust axis, and 
the other is linearized sphericity discussed above. The cut in multiplicity is used partly to 
ensure that the more carefully designed cut indeed perform better. Nevertheless it enhances 
the rate of events with little perturbative activity, as the emitted gluons in general give rise 
to high multiplicities. One advantage of this cut, which makes it worth consideration, is that 
p± observables are left relatively unbiased. 

Fig |l] shows efficiency and purity results for a k±o= 2 GeV for the simulated gluon 
emissions. Most two-jet cuts perform similarly, but it could be noted that the commonly 
used cut in 1 — T not obviously is the most successful. For large rapidities of the emitted 
gluon jet, 1 — T is approximately Ak\/ s and thus not very restrictive in kj_. In the central 
rapidity region, however, 1 — T behaves as k±/y/s. For a study in this region, a cut in 1 — T 
therefore performs very similarly to the other ones considered in Fig |I|. The multiplicity cut 
performs noticeably worse than the others. This is expected, since it is not an cut explicitly 
in a p± variable. 

We end this section by noting that the purity and efficiency measures are not very 
sensitive to the choice of cut. There is thus no optimal cut in all respects, and as we will 
see, the preferred cut instead depends on the observable under investigation. 

3 Observables 

In this section we want to study a set of observables which depend on the properties of 
the hadronization mechanism, but where this dependence is masked by perturbative gluon 
emission. Thus we will not discuss here the flavour composition of the hadrons, as this 



3 




Figure 1: Purity and efficiency for different two-jet cuts. 



appears to be less affected by gluon radiation. (We note, however, that this may not be 
completely true, as there are experimental indications for a larger procution rate of A- 
particles in T decays |fL3 |, and also a larger rate of high-energy ^-particles in gluon jets [ 14] , 
than expected from available fragmentation models.) 

The observables discussed here are all related to transverse momentum. In many mod- 
els p± of neighbouring hadrons are anti-correlated to a larger or lesser degree. This anti- 
correlation is counteracted by jet (or minijet) emission, which gives a positive correlation 
between hadrons close in rapidity, making the analysis of the soft p± generation very non- 
trivial. 
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3.1 Inclusive ^-distribution 

The models discussed in the following sections all have tunable parameters for the soft 
p± generation. The values of these parameters can be adjusted to reproduce the inclusive 
distribution of the hadronic p±. Thus this distribution acts as a constraint and cannot be 
used to separate the different models. For this we have to use more complicated observables, 
e.g. collective variables or correlations. 

3.2 transfer 

One collective variable which can distinguish between different models is the vector sum 
Q_i_ (y) of all hadrons with rapidity less than a given rapidity y 

QX(V) = E P-L<: (3) 

m<y 

which measures the p± transfer over the rapidity y. (This variable has frequently been 
advocated by E. de Wolf |15|]-) F° r a longitudinal phase space model constrained by overall 
momentum conservation, Q± grows with the number of particles as a random walk. Thus in 
this model Q± is larger in the central region and larger in events with high multiplicity. In 
many fragmentation models, however, the p_i_ of a hadron is compensated locally in phase 
space. In such a model, Q± is limited and essentially independent of rapidity and multiplicity. 

An essential point when using the Q± variable is the choice of axis. In an experiment the 
direction of the initial qq pair is not known. If the thrust axis is used to define the transverse 
directions, Q±{y) is kinematically constrained to be equal to zero for y = 0. This definition 
is therefore very unsuitable, if one wants to use the value of Q± in the central region to 
distinguish between different models. 

We have measured Q± w.r.t. different axes in MC generated events, using the Jetset 
default fragmentation scheme, where p± is compensated essentially by neighbouring particles 
in rapidity. This implies that we expect a rapidity plateau in Q±, though decays of unstable 
hadrons and the perturbative activity present in the analysed two-jet events introduce cor- 
rections to this result. As seen in Fig both the thrust axis and the linearized sphericity 
axis by construction create a significant dip in Q± when y = 0. The sphericity axis, on the 
other hand, reproduces the plateau extremely well. This should not be over-interpreted, for 
the reasons mentioned above, but it is clear that the sphericity axis is to prefer in a Q± 
analysis, and it will be used in the results presented below. 

The collective quantity Q± is sensitive to correlations between neighbouring hadrons. It 
is also possible to study directly the correlation in p± between two hadrons separated by a 
rapidity interval 5y 

(PxiVi) ■ P±(V2 = Vi + Sy)) , (4) 
or the correlation between for two separated y-values 

(Q±{yi)-Q±(y2 = yi + 6y)). (5) 

The interpretation of these observables will however be hampered by the systematic effects 
related to the reconstruction of the event axis. In Fig [|we compare (Qj_ ■ Qj_) (5y) for two 
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Figure 2: The properties of (Q\) in MC simulated events depend on the choice of axis. The 
thrust axis, which maximizes longitudinal momentum, by construction gives Q± = when 
y = 0. A similar behaviour is seen for Linearized sphericity. The conventional Sphericity 
axis, minimizing squared transverse momentum, gives the central rapidity plateau expected 
from the chosen string fragmentation model. 
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Figure 3: An explicit measurement of the correlation length Sy is distorted by the choice of 
axis. Two very different p± compensation models are compared when measured w.r.t. the 
original qq axis (left) and the sphericity axis (right). To emphasize the effects of axis choice, 
both the parton cascade and hadron decays have been turned off in the simulations. 



models, one with local p± conservation and one with essentially uncorrelated p±. In the MC 
simulations, we have taken the liberty to turn off both the parton cascade and the hadron 
decays, in order to more clearly illustrate the effects of axis choice. The plots show the results 
obtained w.r.t. the original qq axis and the sphericity axis, respectively, and it is clear that 
the behavior of (Qj_ ■ Q_l) (Sy) is drastically changed. (Equally drastic changes occur if using 
the thrust axis or linearized sphericity axis.) Due to these effects of axis reconstruction in 
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e + e~, a variation of 8y does not seem to provide much more information than is available 
from the observable (Q\). This orgument does not imply that the Sy correlation cannot 
be an interesting observable in deep inelastic scattering and hadron-hadron collisions, where 
the event axis is better known. 



3.3 Screwiness 

A particular correlation between rapidity and azimuthal angle is determined by a variable 
called screwiness and introduced in ref ||. Screwiness is defined by 



S(u>) = 



(N(y cut )) 



exp(i(Luyi - <&)) 

|»i[<lteut 



(6) 



where yi and (pi is the rapidity and azimuthal angle for particle i, and the sum runs over all 
central rapidity particles, specified by some y cn f For large u the sum describes a random 
walk along unit vectors in the complex plane, and the absolute square of the sum becomes 
the average multiplicity in the considered interval, (N(y cut )). We have normalized S (u>) by 
this quantity, which implies that it is 1 for large uj. 

The screwiness measure is constructed to give a signal for a helix-like correlation between 
y and <fi. A motivation why such a correlation might be expected is presented in ref || and 
shortly described in the following section. 



4 Models 

The most studied hadronization models are based on string dynamics or cluster fragmen- 
tation. In the cluster model implemented in the HERWIG MC ||T^[ the gluons split into qq 
pairs at the end of the cascade. A quark is then combined with an adjacent antiquark to 
form a colourless cluster, which normally decays to two (possibly unstable) hadrons. This 
approach becomes less realistic and must be modified in events with unusually low gluon 
radiation. In such events some clusters can get very large masses, and in the Herwig MC 
these clusters decay in a stringlike fashion. Consequently, if we select events with low gluon 
activity we cut away those events for which Herwig is meant to work best. For this reason 
we feel that it would not be meaningful to compare such an event sample with Herwig 
simulations, and therefore we will instead in the following study some different versions of 
string hadronization, and for comparison also a model with a p± behaviour similar to a 
longitudinal phase space. 



4.1 String fragmentation 

In the Lund string model the confining colour field is assumed to behave as a massless 
relativistic string starting at a quark, passing the gluons in accordance with their order in 
colour, and ending at an antiquark. Thus the gluons act as transverse excitations or kinks 
on the string. The string breaks by the production of qq pairs (or occasionally (qq)(qq) 
systems) which combine to hadrons. 
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For the longitudinal momentum distribution it is assumed that any breakup divides the 
string in two causally disconnected pieces, which decay further independently of each other. 
Together with the assumption of a central rapidity plateau this implies a unique result, the 
symmetric Lund fragmentation model |17], [UJ. The result can be generated in an iterative 
way with a splitting function of the form 

/Woc iI^ eX p ( _ 6M |/z), (7) 

where Mj_ is the transverse mass of the generated hadron and z is its fraction of jet energy 
plus longitudinal momentum. 

The iterative procedure specifies an ordering of hadrons called rank. The original quark 
qo combines with an antiquark qi from the qiq x breakup into the rank one meson, qi then 
forms the rank two meson together with the antiquark q~2 from the next breakup, etc. Though 
the order in rank is not an observable, it is strongly correlated to the order in rapidity. 

The string breakup by qq creation has been treated as a tunneling process [18|, in 
which the transverse momenta should be described by a Gaussian with a width determined 
by the string tension k: 

-^-j oc exp(-m 2 ± /a 2 ). (8) 

Here m 2 ^ = m 2 + k\ is the squared transverse mass of the quark or the antiquark, and from 
the tunneling mechanism the width should be given by a 2 = k/tt. 

An essential property of the string fragmentation is infrared stability, meaning that a 
soft or collinear gluon has only little influence on the string motion and therefore also on 
the final hadronic state []• This feature makes the result less sensitive to the cutoff for the 
partonic cascade, and a change in the cutoff can be well compensated by a change in the 
hadronization parameters when looking at event shapes and inclusive observables for the 
total event sample ||. Thus a larger cutoff is compensated by a larger width a for the 
hadronization k±. Some differences may remain, however, as a gluon tends to give extra 
k± in the same direction to neighbouring hadrons, while hadronization-p^ tends to be in 
opposite directions, depending on the details in the fragmentation model. 

In section we discussed an observable called Q^, which is sensitive to the properties 



of p± compensation. We will continue this section by describing different hadronization-p^ 
models, which give slightly different predictions on Q± and the other observables discussed 
in section |3]. In section |5] we examine how well these differences can be distinguished in a 
two-jet event sample. 



4.2 Jetset default 

In the Jetset default version of the string fragmentation model the qq breakups of the 
string are uncorrelated both in transverse momentum and in flavour. As the quark and 
the antiquark in a pair have balancing k± the transverse momentum of one hadron is fully 

3 To treat situations where a hadron obtains energy and momentum from several small string pieces is 
however not uniquely specified from the string dynamics. In the models studied below we will use the method 
developed by Sjostrand (2(J and implemented in the Jetset MC. 
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compensated by its two neighbours in rank. This implies that p± is conserved very locally 
in rapidity. 

From Eq (J|) we expect that the width of the /^-distribution should be given by y 1 n/n « 
250 MeV, but in applications it is treated as a tunable parameter, correlated to the cut-off for 
the parton cascade. Fits to data give a value around 360 to 400 MeV. Eq (||) also introduces 
an overall suppression of heavy quarks. This implies that the rate of c is completely negligible 
(the c/u ratio is ~ 10~ n ). To get the rate of s-quarks we would have to know its effective 
mass. The observed s/u ratio around 1/3 corresponds to the very reasonable value 250 MeV. 

Studies of photons at LEP1 indicate a larger tail for high-p^ tt°s [ 2X[| . The origin of this 



tail is not clear, but in Jetset default a small component with higher p± is included (k± is 
multiplied by a factor of 2 for 1% of the vertices). 

4.3 Partial p ± Compensation 

The production rate for rj and rf cannot be reproduced if neighbouring string breakups are 
uncorrelated in flavour. Also studies of two-particle correlations observed in hadron-hadron 
collisions are not well described if the in different breakups are totally uncorrelated [^2[ . 
As discussed above, this would imply that the p± of a hadron is fully compensated by its 
two neighbours in rank. In ref a model is presented in which there are correlations in 
flavour and k± with some finite correlation length. The model implies a suppression of 77 
and rj' production, and also that only a fraction 7 < 1 of the p± of a hadron is compensated 
by the neighbouring hadrons in rank. Thus p± is compensated within a finite correlation 
length, but this range is larger than in Jetset default. 



In it is shown that this p± compensation assumption can be embedded in an iterative 



fragmentation scheme. The inclusive p± distribution is then a Gaussian, 

dP cx exp(-— pl)d 2 p ± , (9) 
27 

where a is a tunable parameter. Furthermore, the sum in for the first i particles in rank, 
Q_i_i, saturates for large i towards the distribution 

dP oc exp(-aQi)d 2 Q ± . (10) 

We note that for fixed (p 2 ± ) for the hadrons, the collective variable (Q 2 ^) is smaller for larger 
7. This is just the previously stated behaviour for the closely related observable (Q 2 ±(y)), 
which is small if p± is conserved locally in rapidity. 

The Jetset default algorithm is retrieved by setting 7 = 1 and a = 1/cr 2 , where a is 



the width in Eq (|8|). In J23|, 7 is instead assumed to depend on the hadron mass M, and is 
parametrized as 

7 = 7—, 11 

' 1 + M /M' v ' 

where M is set to 0.9 GeV. 
4.4 Uncorrelated pj_ 

To test whether the data really show a local p± conservation, we also want to study a 
situation where the hadronic transverse momenta are uncorrelated apart from overall p± 
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conservation. This can be achieved by pushing M to infinity in Eq flTTD, which corresponds 
to 7 being small, but keeping a/7 finite. The MC from is not designed to cope with very 
large M -values. When hadrons are peeled off from both ends of the string usually a large 
p± is built up for a remaining string piece in the centre. Such events have to be discarded, 
keeping only those events where by chance this string piece has limited p±. In this paper we 
study only central particles in the accepted events, and thus for our purposes we obtain a 
fair estimate of the result if in the MC we peel off hadrons from one end only, pushing the 
junction to the other and outside the region of study. 



4.5 The UCLA model 

A partial p± compensation is also assumed in the UCLA model ||. In this model it is assumed 
that not only the longitudinal momentum distribution but also the flavour composition and 
the distribution in transverse momentum are determined by the symmetric fragmentation 
function in Eq (0). Thus in this model the quark degrees of freedom do not appear explicitly, 
but only the momenta and flavour of the hadrons. A direct application of Eq (|7]) will however 
give a too wide ^-distribution. In the UCLA approach this problem is solved by assuming 
a similar kind of finite p± correlations as in the model in section ^]3|, with 7 = 1/2. Thus 
as far as transverse momentum goes we expect rather similar results for the UCLA model 



and the model in [P3|1 , and for this reason we content ourselves to study the latter in our 
investigatons. 



4.6 Helix String 

At the end of a perturbative parton cascade the virtualities become relatively small and 
the running coupling therefore relatively large. This implies that interference and coherence 
effects are expected to be important. In ref arguments are presented for correlations 
between rapidity and azimuthal angle for the emitted gluons, which correspond to a helix- 
like field at the end of the cascade. This field configuration has also some similarities with 
a string coupled to a gauge field or a massive relativistic string [^4| . 

A modification to the Lund fragmentation scheme reflecting the helix-like field is pre- 
sented in [[J. The rapidity distance Ay between two string breakups is correlated with the 
azimuthal angle A0 between the transverse direction of the two produced qq pairs. The pitch 
of the helix is given by an unknown parameter Ay/A<p = r, allowing for some Gaussian fluc- 



tuations around this value. This implies that the observable S(u) defined in section |3.3| gets 
a peak at uj « 1/r. Expected values of r are around 0.3 to 0.5. The string energy per unit 
rapidity available to generate p± is given by a parameter m, which implies that a typical 
magnitude of hadron p± is given by mr. 



4.7 Time Dependent p± 

In order to investigate the sensitivity of the observables we will also study a model with 
a correlation between transverse momentum and multiplicity which might not be totally 
inconceivable. The results in Eq (^) with a constant width corresponds to a situation with 
a static flux tube or stringlike colour field. The estimated width of this flux tube is of the 
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order of 1 fm. One can imagine that at very early times the flux is more concentrated. This 
might correspond to a larger effective string tension, and also from the uncertainty relation 
one might then imagine that larger transverse momenta are produced in early breakups. If 
the transverse dimension of the flux tube equals the proper time r, then we could expect 
fluctuations in k± of the order 2/r. Expressed in the variable T = k 2 t 2 this corresponds 
to a width a 2 ~ Ak 2 /T. This becomes very large when r — > 0, and from a string picture 
we should expect a saturation. To produce large kj_ sufficient energy must be stored in 



the string, which therefore must have a minimum length. As discussed in |TjJ we expect 
a suppression when k± > kt, and thus a saturation when 2/r ~ nr or T ~ 2k. These 
arguments would give the following width for the fc_|_-generation 

"M^*-^ rt% (12) 

I %ax = - + 2k, r > 2k 

Actually this expression does reproduce the inclusive p± spectrum with hardly any tuning. 
In our calculations we allowed for an overall normalization factor on a, 

a -> xcx. (13) 

For the standard perturbative cutoff we get \ = 1.08. 

A special consequence of the time dependent kj_ width is a strong correlation between 
transverse momentum and multiplicity. An early breakup splits the original system in two 
pieces with relatively small invariant masses. If these masses are called Mi and M 2 we have 
approximately 

M\M\ « Ts, (14) 

where a/s is the total mass of the system. It is also easy to show that the ensuing depletion 
in multiplicity appears in a rapidity of size Ay ~ ln((r) /T) around the breakup, where (T) 
denotes the average value of T. 



5 Results 

5.1 Hadron p± spectra 

We have used average p± and average multiplicity to impose constraints on the free param- 
eters of the considered models. For simplicity, we have left the Ariadne+Jetset default 
algorithm unchanged and tuned the other models to reproduce its values. Differences are 
then searched for in other observables, which will be discussed below. 

All models depend on the k± cut of the cascade and the a and b parameters of the symmetric 
fragmentation function Eq (0). To compare different hadronization models, we are interested 
in using the same description of the perturbative phase for all models, and we have therefore 
tuned all models keeping the cascade cut-off k± cut constant to the Ariadne default value. 
As fcj_ C ut is also an hadronization parameter, we will in some examples compare tunes of the 
same model with different values of fcj_ cu t- 

The a and b parameters influence the multiplicity distribution. For each setup of other 
parameters, there is in general a narrow but long band of values in the a and b parameter 
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model 


parameters 










(N) 


(P±) 




k± c ut 


a 


b 


others 












(GeV) 




(GeV -2 ) 












Default 








a k± (GeV) 












0.6 


0.23 


0.34 


0.405 






34.9 


0.60 




1.5 


0.5 


0.4 


0.42 






35.1 


0.59 


Partial p± 








75 (GeV) 


M (GeV) 








Compens. 


0.6 


0.23 


0.34 


0.67 


0.9 




34.9 


0.60 


Uncorr. p± 


0.6 


0.23 


0.34 


1.52 


10 




34.8 


0.59 


Helix 








T 


a T 


m(GeV) 








0.6 


0.23 


0.37 


0.3 


0.2 


1.1 


35.0 


0.59 




0.6 


0.23 


0.39 


0.5 


0.3 


0.8 


34.8 


0.60 




0.6 


0.23 


0.4 


0.7 


0.35 


0.68 


34.6 


0.60 


Time 








X 










Dep. p± 


0.6 


0.23 


0.34 


1.08 






34.7 


0.60 




1.5 


0.5 


0.42 


1.18 






34.7 


0.60 



Table V. The tuned parameter values of the considered models. The tunes are performed 
w.r.t. average multiplicity (N) , with 7r° stable, and average (p±) within a rapidity interval 
\y\ < Vent = 3. (p±) is measured w.r.t. thrust axis and defined as ncv{1 ^ iycut)) Ecv Y,\ y \<y cut \P±\- 
The top line of the models is the Ariadne+Jetset default tune which has not been changed. 
The other models are tuned to give similar results. The models and their parameters are 
discussed in section [|. 



space which reproduces the average multiplicity. The specific a and b values in the tunes 
are chosen to give a fair agreement also on multiplicity dispersion. The result is presented 
in Table |I[ 



p± in low multiplicity events 

The different models have been tuned to give similar p± distributions in the full event sample. 
After a two-jet cut, the p± distribution may differ, as the importance of the perturbative 
phase is then reduced. We will here investigate that possibility. Among the two-jet cuts 
discussed in section |2], the rather blunt cut in multiplicity is then of interest, as it introduces 
little bias to the p± spectrum. 



In section |4~7| , we discussed a model where the proper time of a qq vertex in the string 
influences the width of the k± distribution for the qq pair. This model introduces correlations 
between average p± and multiplicity. 

Gluon emissions in general contribute both to p± and multiplicity. For soft gluons the 
effects are smaller, but if their contribution to p± dominates, it affects the average p± in 
low-multipicity events. Thus it is of interest to compare correlations in (p±) and multiplicity 
for the same fragmentation model tuned with different cascade cutoffs. 

In Fig [| we compare the time dependent p± model discussed in section |4.7| and Jetset 
default fragmentation by studying average p± of particles in the central rapidity region 
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Figure 4: Average p± w.r.t. thrust axis for different charged multiplicities. Statistical errors 
are well within symbol sizes. The correlation between p± and multiplicity in the time depen- 
dent p± model is seen as a higher prediction for average p± in low multiplicity events. For 
both models (p±) is higher in low multiplicity events when the cascade cutoff k± cnt is lowered. 
This indicates that soft gluons in these models give a larger contribution to hadron p± than 
to hadron multiplicities. 

\y\ < 2, as a function of the charged multiplicity N c ^ in that region. As expected, average p± 
at low multiplicities is larger in the time dependent p± model. p± and y are defined w.r.t. 
the thrust axis, and no two-jet cut in a p± measure is imposed on the event sample. 

We note that (p±) at low multiplicities is higher for both models if the cascade cutoff is 
reduced from 1.5 GeV to 0.6 GeV. In the dipole cascade combined with string fragmenation, 
the soft gluons in this range of k± thus contributes to hadron p± more than to hadron 
multiplicities. 

We have also checked the prediction for the time dependent p± model if the expected 
saturation of the qq k± width at very early breakups is neglected. We then get a significantly 
higher prediction for (p±). When iV c h(|y| < 2) = 2, it lies between 0.9 and 1.2 GeV, 
depending on the k± cnt for the cascade. The correlation in p± and multiplicity is thus 
and observable which can be used to confront this kind of unexpected p± generation with 
experimental data. 



5.2 p± transfer 



The observable Q±(y) discussed in section [372] measures the p± transfer over the rapidity 
y and is sensitive to correlations in hadron p±. We have investigated Q± for different as- 
sumptions about hadron p± conservation. In the default algorithm, p± is compensated by 
the neighbours in rank, which implies a compensation local in rapidity. In the partial p± 
compensation model, a fraction 7 < 1 of p± is compensated by closest neighbours. This 
implies a correlation length for p± of order I/7 which is finite, but larger than for Jetset 
default. We also examine a model where the hadron p± are uncorrelated, apart from global 
momentum conservation. 
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(Qi)(GeV 2 ) 


default 


partial p± 


uncorr. p± 








compens. 




k± cu t (GeV) 


1.5 


0.6 


0.6 


0.6 


No cascade 




0.60 


0.71 


0.74 


With cascade, all events 


12.2 


12.3 


12.4 


12.6 


two-jet events 


0.86 


0.94 


1.02 


1.09 



Table 2: Average Q\, measured w.r.t. the sphericity axis in a central rapidity range \y\ < 2. 
Results for different assumptions about hadron p± correlation lengts. The differences in the 
predictions are shadowed by the perturbative cascade, but can be restored by a two-jet cut. 
Here a cut in sphericity, ln(S') < —4.5, has been used. 



In two-jet events, the hadrons can also aquire p± from soft gluons. The compensation of 
this p± is in general not identical to the one assumed in the fragmentation, and depends on 
the recoil treatment of the cascade formalism. Thus the assumed p± conservation length may 
depend on the treatment of recoils in the cascade, and on the cutoff scale, which determines 
to what extent hadron p± originates from the cascade or the hadronization. 

The top row in Table |] shows results without cascade for the three different p± correlation 
assumptions. Q± clearly grows for less local p± compensation, with a relative difference of 
25% between complete p± compensation by neighbours (Jetset default) and uncorrelated 
p±. Adding a cascade allmost wipes out the difference, but after a two-jet cut it is restored to 
a satisfactory degree. The relative difference between the two extremes with a common cutoff 
k± C ut= 0.6 GeV (second and fourth column) is then about 16%. Also partial p± correlation 
as assumed in p3[ , similar to the UCLA p± generation model, gives results after a proper 
two-jet cut which differ significantly from the results of the Jetset default assumption. 

A comparison of the first and second column illustrates the difference in p± compen- 
sation in the cascade and in the fragmentation. The same p± model, tuned with different 
cascade cut-offs, give significantly different results on Q±. Soft gluon emission in the dipole 
cascade introduces longer p± correlation lengths than given by default hadronization p±, 
which implies that Q± increases when adding soft gluons in the range 1.5 to 0.6 GeV to the 
cascade. 



5.3 Screwiness 

In ||, MC simulations of qq string fragmentation without cascade show a clear signal in 
S(u) in the helix-like fragmentation model for r ~ 0.5 or larger. To check the influence of 
relatively soft gluons on the screwiness measure, we have modified Jetset to generate p± 
according to the helix model. The p± generation is then combined with the Jetset default 
recipe for passing gluons f2~0fl . This implies that the azimuthal directions on different sides of 



a gluon kink is not completely uncorrelated, though measured w.r.t. two different directions. 
The p± component being transverse to both considered directions is preserved when passing 
a gluon kink. 

One could also consider more direct effects, assuming the azimuthal angle of vertices on 
different string segments to be explicitly uncorrelated, and also allowing the helix to have 



14 




Figure 5: Screwiness without gluons, and with low-k± gluons (k± in the emissions is re- 
stricted between 1 and 2 GeV). The central rapidity region \y\ < 2 is considered. With the 
assumptions made here about effects from perturbative gluons, the signal survives. 



different orientation on different string segments. The investigation here however focus on 
the minimal effects of gluons bending the string. 

Before looking at screwiness after a full cascade and a two-jet cut, we examine the effects 
of soft gluons by running a cascade with k± for the emissions constrained between 1 and 2 
GeV. As seen in Fig |5|, the observable S(u) is fairly unaffected by these rather soft gluons. 

A cascade introduces long-range azimuthal correlations which have large effects on S(u). 
In a three-jet event the thrust direction tends to be aligned with the hardest jet, which 
implies that the particles of this jet are rather randomly distributed in azimuthal angle (f>. 
The softest jet is seen as a set of particles well collimated both in (f) and y. The particles 
of the remaining jet are more spread out towards large rapidities and have and azimuthal 
direction opposite to the softest jet. In this typical three-jet configuration, particles of the 
hardest jet merely introduces noise to S(u), while the two softer jets tend to enhance S(u>) 
when uj ~ 7i/y cnt , where y cut is the largest considered rapidity in the analysis. Thus we get a 
signal which is independent of fragmentation mechanism, but is sensitive to the magnitude 
of the considered central rapidity range. This behaviour is seen in the upper plots of Fig || 

After event cuts (lower plots in Fig ^), the result is very similar to those with low-A;_i_ 
gluons in Fig [5] and give a clear signal for helix-like string fragmentation if r ~ 0.5 or larger. 
In 0, different ways to enhance the signal also for lower values of r were investigated. The 
analysis was performed on qq strings without cascade, and it was shown that a lower limit on 
average p± could be used to distinguish default string fragmentation and helix fragmentation 
with t down to 0.3. A lower cut in p± is slightly in conflict with a two-jet cut, but it is 
possible that a similar approach can be used to enhance the signal in two-jet events obtained 
from real data. 
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Figure 6: Screwiness after gluon cascade, with and without a two-jet cut. Without any 
two- jet cuts (upper row), a signal reflecting long-range azimuthal correlations in three-jet 
structures is seen, independent of fragmentation mechanism. After a two-jet event cut, this 
signal is removed, and the S(u) measure is sensitive to helix-like fragmentation. The results 
are similar to the ones with only low-k± gluons, and only moderately reduced. 

6 Conclusions 



High energy reactions like e + e~ — > hadrons are usually described in terms of two phases, an 
initial perturbative phase, formulated as a parton cascade, followed by a soft hadronization 
phase described by a phenomenological model. It is however not only the description of the 
soft phase which depends on model assumptions. This is also the case for the perturbative 
phase towards the end of the cascade, where a s is large and some kind of cut-off is neces- 
sary. The transition region may also exhibit interesting interference or coherence effects, one 
example being the helix-like field configuration proposed in |J. 

In many cases the effects of gluon emission overshadow the features of the hadroniza- 
tion phase. To investigate the hadronization mechanism it is therefore interesting to study 
whether suitable events cuts can help separating the two phases and give observables sensitive 
mainly to one of them. 
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We have examined different strong two-jet cuts to suppress the perturbative activity. It 
is possible to define efficiency and purity measures for the different cuts. MC simulations 
show, however, that no cut can be called optimal. Instead the preferred cut depends on the 
observable under investigation. 

As the flavour composition of the final state in most models is less sensitive to the parton 
cascade, we concentrate here on observables related to transverse momentum. Different 
assumptions can be made concerning recoil effects and the locality of p± conservation. In 
most models it is possible to tune the inclusive p± spectrum, and to distinguish the models 
it is then necessary to study different forms of correlations. The vectorial sum of for all 
particles with rapidities smaller than some value y, here called Q^, is an observable which 
is sensitive to the p± correlation length. If p± is locally conserved, average Q± is similar to 
average p±. On the other hand, if only global constraints restrict p±, the collective measure 
Q± gets significantly larger than p±. 

After a two-jet cut, we find the Q± measure to be sensitive to assumptions about p± 
compensation in the hadronization phase. It is also sensitive to the value of the cascade 
cut-off, since the recoil treatment in the cascade in general give longer p± correlation lengths 
than assumed in the fragmentation. 

Another observable which is shown to be sensitive to the cascade cut-off is average p± 
in low multiplicity events. Though multiplicity in general is not considered a two-jet cut, it 
suppresses high-fc^ perturbative emissions, which typically give rise to large multiplicities. 
Gluon emissions with moderate k± however contribute more to hadronic p± than to hadronic 
multiplicity, which implies that the correlation between p± and multiplicity is sensitive to 
the cascade cut-off. 

In || it is argued that coherence effects in the transition region may give rise to a helix- 
like colour field, with a correlation beteen rapidity and azimuthal angle. Perturbative gluon 
emissions mask this correlation to a large extent, but we have found that after a strong 
two-jet cut, a signal in the proposed observable "screwiness" is present for some possible 
helix configurations. 

The very high statistics now available from experiments at the Z° resonance implies 
that strong two-jet cuts are possible on the data without loosing the statistical significance. 
Our results show that such cuts can indeed in many cases discriminate between different 
hadronization models and thus be a tool for a more detailed study of the hadronization 
mechanism. 

References 

[1] D. Amati, R. Petronzio, G. Veneziano, Nucl. Phys. B140 (1978) 54; 
ibid. B146 (1978) 29 

[2] Ya.I. Asimov, Yu.L. Dokshitser, V.A. Khoze, S.I. Troyan, Z. Phys. C27 (1985) 65 

[3] B.R. Webber, Nucl. Phys. B238 (1984) 492 

[4] CD. Buchanan, S.B. Chun, Phys. Rev. Lett. 59 (1987) 1997; 
CD. Buchanan, S.B. Chun, Phys. Rep. 292 (1998) 239 



17 



[5] B. Andersson, G. Gustafson, A. Nilsson, C. Sjogren, Z. Phys. C49 (1991) 79 

[6] B. Andersson, G. Gustafson, J. Hakkinen, M. Rigner, P. Sutton, JEEP 09 (1998) 014 

[7] L. Lonnblad, Comp. Phys. Comm. 71 (1992) 15 

[8] M. Bengtsson, T. Sjostrand, Comp. Phys. Comm. 39 (1986) 347; 
T. Sj6strand, Comp. Phys. Comm. 82 (1994) 74 

[9] G. Gustafson, Phys. Lett. B175 (1986) 453; 

G. Gustafson, U. Pettersson, Nucl. Phys. B306 (1988) 746; 

B. Andersson, G. Gustafson, L. Lonnblad, Nucl. Phys. B339 (1990) 393 

[10] B. Andersson, G. Gustafson, G. Ingelman, T. Sjostrand, Phys. Rep. 97 (1983) 31 

[11] L. LSnnblad, Z. Phys. C58 (1993) 471 

[12] G. Parisi, Phys. Lett. B74 (1978) 65; 

J.F. Donoghue, F.E. Low, S.Y. Pi, Phys. Rev. D20 (1979) 2759 

[13] W. Hofmann, Ann. Rev. Nucl. Part. Sci. 38 (1988) 279 

[14] L3 Collaboration, M. Acciarri et al, Phys. Lett. B371 (1996) 126 

[15] E.A. de Wolf, private communications. 

[16] G. Marchesini, B.R. Webber, M.H. Seymour, G. Abbiendi, L. Stanco, I.G. Knowles, 
Comp. Phys. Comm. 67 (1992) 465 

[17] B. Andersson, G. Gustafson, B. Soderberg, Z. Phys. C20 (1983) 317 

[18] H. Bohr, H.B. Nielsen, NBI-EE-78-3 (1978) 

[19] B. Andersson, G. Gustafson, T. SjSstrand, Z. Phys. C6 (1980) 235 
[20] T. SjSstrand, Nucl. Phys. B248 (1984) 469 

[21] ALEPH Collaboration, D. Buskulic et al., Z. Phys. C57 (1993) 17; 
L3 Collaboration, O. Adriani et al., Phys. Lett. B292 (1992) 472 

[22] E.A. de Wolf, Contrib. to Proc. XXE Symp. on Multiparticle Dynamics, Santiago de 
Compostela 1992, ed. A. Pajares, World Scientific 

[23] B. Andersson, G. Gustafson, J. Samuelsson, Z. Phys. C64 (1994) 653 

[24] N.K. Nielsen, Nucl. Phys. B167 (1980) 249 



18 



